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Abstract 

We describe the effects of phase coherence on transport and thermodynamic prop- 
erties of a disordered conducting network. In analogy with weak-localization cor- 
rection, we calculate the phase coherence contribution to the magnetic response of 
mesoscopic metallic isolated networks. It is related to the return probability for a 
diffusive particle on the corresponding network. By solving the diffusion equation 
on various types of networks, including a ring with arms, an infinite square net- 
work or a chain of connected rings, we deduce the magnetic response. As it is the 
case for transport properties -weak-localization corrections or universal conductance 
fluctuations- the magnetic response can be written in term of a single function S 
called spectral function which is related to the spatial average of the return prob- 
ability on the network. We have found that the magnetization of an ensemble of 
connected rings is of the same order of magnitude as if the rings were disconnected. 
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I. INTRODUCTION 



The understanding of persistent currents in mesoscopic rings is still an open prob- 
lem since no quantitative agreement has yet been reached between theoretical de- 



velopments and experimental results 


Buttiker et. al. 1983, Levy et al. 1990 


, Chan- 


|drasekhar et al. 1991|JMailly et al. 1993, 


Keulet et al. 1995 


: Mohanty et al. 1996, 


LVlontam- 



[baux 1996|| . The only quantitative calculations of persistent currents relevant to the 
experimental situations rely on analytical methods in which disorder and interactions 
are treated perturbatively [|Cheung et al. 1989| , |Ambegaokar et al. 1990|J3chmid 1995| , |Von 



Oppen et al. 1991|,|Altshuler et al. 1991 


Akkcrmans et al. 1991 


,|Oh et al. 1991 


Arga- 



man et al. 1993,Argaman et al. 1993,CTontambaux 199b'|. The disorder is described in 



the usual framework of perturbation theory and interactions are treated in a Hartree-Fock 
approximation. 

In ring experiments, the important parameters are certainly the strength of disorder 
and the strength of interactions. Both parameters are difficult to vary experimentally in a 
controlled way. In order to go beyond this problem, we have recently proposed to extend 
the analytical method to describe new geometries which generalize the ring topology, and 
which allow the variation of some geometric parameters in order to change the amplitude 
of the persistent current | |Pascaud et al. 1997| |. In the experiment done on a single Ga- 
As ring, for example, the ring was connected to leads ||Mailly et al. 1993H and one may 
wonder how the current is changed when the lengths of the leads increase. On the other 
hand, experiments on many rings have been performed on an ensemble of disconnected rings 
Levy et al. 1990| , [Reulet et al. 1995|| . One may wonder whether the same effect persists 



when the rings are connected, i.e. what is the magnetization of a lattice, resulting from 
phase coherence. 

Such a generalization, from the ring geometry to a network, had already been considered 
in the past to describe the weak-localization correction to the conductivity. This correction 
is related to the return probability, solution of a diffusion equation [[Khmclnitskii 1984 ] . This 
return probability has two components, a purely classical one and an interference term which 
results from interferences between pairs of time-reversed trajectories. In the diagrammatic 
picture, they are related to the diffuson and Cooperon diagrams. The interference term, 
p 7 (r, t',uj), is field dependent and is solution of the diffusion equation: 



[ 7 _ _ D(V + ^) 2 K(r, r\ w) = *(r 



(1) 



where D is the diffusion coefficient. The scattering rate 7 = D / 'lA describes the breaking 
of phase coherence. is the phase coherence length. 7 will be compared to 1/td = D/L 2 
where To is the diffusion time, typical time to diffuse through the system of size L. This time 
is the inverse of the Thouless energy. The classical probability obeys the same equation as ([]]) 
with A = 0. By solving the diffusion equation (JJ) on various lattices, the weak-localization 
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correction was calculated by Dougot and Rammal [ Poucot et al. 1986| . The persistent cur- 
rent of a ring, or in a more general geometry the mesoscopic magnetization, have a different 
nature. They are thermodynamic and not transport quantities. Transport properties exhibit 
weak-localization effects because they are response functions related to the propagation of 
two states: the conductivity is written in terms of a product of two Green functions whose 
disorder average is related to the return probability. A priori, average thermodynamic quan- 
tities do not exhibit weak-localization corrections since they are integrals of the density of 
states, i.e. they are written in terms of a single Green function. This is why the average 
persistent current (/) of non interacting electrons vanishes for an ensemble of rings. How- 
ever, the typical current, I typ = J (I 2 ) involves the product of two propagators. Similarly, 
the average Hartree-Fock correction to the average current does not vanish because the in- 
teraction term also contains two propagators. As a result, these quantities are expected to 
be related to weak-localization effects and to the return probability. 

This paper is organized as follows: in the next section, we recall how to solve the diffusion 
equation on a lattice and we calculate the weak-localization correction. In section III, the 
persistent current and the mesoscopic magnetization are related to the return probability 
and are calculated for a few geometries in section IV. In the last section, we show how the 
universal conductance fluctuations are related to the return probability and we calculate 
them for the simple geometry of the single wire, using this simple approach. 



II. DIFFUSION ON A LATTICE - WEAK-LOCALIZATION CORRECTION 

Let us first recall that the weak-localization correction to the conductance of a connected 
mesocopic sample can be written in term of the interference part of the return probability 
Khmelnitskii 1984jjDoucot et al. 1986U : 



Acx(r) = -2s(e 2 /h)DC 7 (r, r) (2) 

s is the spin degeneracy. The Cooperon C 7 (r, r, H) is the time integrated field- dependent re- 
turn probability: C 7 (r, r, H) = / °°p 7 (r, r, t, H)dt = p 7 (r, r, uj = 0, H). Dougot and Rammal 



Doucot et al. 1986] have calculated the weak-localization corrections on various networks. 
Considering networks made of quasi- ID wires, so that the diffusion can be described as 
one-dimensional, the Cooperon C 7 (r, r') obeys the one-dimensional diffusion equation 

[ 7 - D(V + -^-) 2 }C,(r, r') = 6{r - r') (3) 

with the continuity equations written for every node a (including the starting point r' that 
can be considered as an additional node in the lattice) ||Doucot et al. 198"6 
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r, r' are linear coordinates on the network. The sum is taken over all links relating the 
node a to its neighboring nodes (3. Integration of the differential equation (|3|) with the 
boundary conditions (f|) leads to the so-called network equations which relate C 7 (a,r') to 
the neighboring C 7 (/3,r'). 



Ecoth(^)C(a,r')-E 



sinh(/ a/3 /L (/) ) DS 



5„ 



(5) 



l a p is the length of the link (a(3) and 7^ = (47r/0 o ) fa Atfl is the circulation of the vector 
potential along this link. Solving this set of linear equations and performing the spatial 



integration of C„ 



(r' r') 



give access to the weak-localization correction. 



III. PERSISTENT CURRENTS AND MESOSCOPIC MAGNETIZATION 



The magnetization M(H) is the derivative of the free energy F with respect to the 
magnetic field H: M — — 4C. Introducing the field dependent DOS (for one spin direction), 
p(e, H), the magnetization is written at zero temperature as (taking the spin into account): 



M = -2 



_d_ 
~dH 



ep(e, H)de 



(6) 



The origin of energies is taken at the Fermi energy. The average magnetization is thus 
related to the field dependence of the average density of states (p(e, H)). In a bulk system, 
this leads to the Landau diamagnetism. In the thin ring geometry where no field penetrates 
the conductor, the DOS only depends on the Aharonov-Bohm flux. Its average is flux 
independent because the flux modifies only the phase factors of the propagator which cancel 
in average. During the last years, there has been a large amount of work to explain the origin 
of a non-zero average current, as observed experimentally. The constraint that the number of 
particles is fixed in each ring leads to an additional contribution to the average current often 
refered as the " canonical" current |[Bouchiat et al. 1989| , [Imry 1991|| . Its amplitude is however 
very small ||Schmid 1995| , |Vo~n Oppen et al. 199T] , [Altshuler et al. 1991| , |Akkermans et al. 1991 



As we will see later, a larger contribution comes from the effect of interactions. 

(M 2 ) — (M) 2 . From 



We first calculate the typical magnetization M typ , defined as M 2 yp 



eq.(P), it can be written as: 



Mi 



d d r u /" u 



tup 



(7) 



where K is the correlation function of the DOS: K(e — e', H, H') = (p(e, H)p(e', H')) — pjj. 



K(e) has been calculated by Altshuler and Shklovskii ||Altshuler et al. 1986|| for a bulk sys- 
tem and later in the presence of a magnetic flux ||Schmid 1995|,[Vo"n Oppen et al. 1991|,|Alt- 



shuler et al. 1991] , [Akkermans et al. 1991|| for the ring geometry. A very useful semiclassical 
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picture has been presented by Argaman et al, which relates the form factor K(t), the Fourier 
transform of K(e), to the integrated return probability P(t) = J p(r,r,t)dr for a diffusive 
particle ||Argaman et al. 1993|| : 



K(t) = tP(t)/(4n 2 ) 



(8) 



In this picture, like the return probability, the form factor is the sum of a classical and 
an interference term: 



K(t,H,H') 



47T 



t r „ , H — H . „ . H -\- H 
[Py{t,—Z—)+P7{t,—Z—)] 



Fourier transforming K(e — e') and using the identity J °° edee ltt 
straightforwardly^ 



Ml P {H) 



P,(t,H)\* 



H 



dt 



(9) 

•1/t 2 , one obtains 
(10) 



2tt 2 Jo t 3 
where P*{t,H)\f? = d 2 P 7 /dH 2 \ H - d 2 P^/dH 2 \ . 

We now turn to the average magnetization. It has been found that the Hartree 



Fock(HF) correction to the energy leads to an average persistent current in a ring [ Am 



[begaokar et al. 1990| , t3chmid 1995] . The average magnetization can be written as ||Arga 

U d 



man et al. 1993|,[Montambaux 1996 



(M ee ) 



dEr 
K dH 1 



4 dH 



(n (r))dr 



(11) 



where the the screened Coulomb interaction is U = 4ire 2 /q 2 ^ F [ bicker n 1991 ] . qTF is the 
Thomas-Fermi vector. The magnetization can be rewritten in terms of the two-point cor- 
relation function of the local density of states and then in term of the return probability 
Montambaux 1996||: 



Upo_d_ r f 7 M) 

7T dH Jo t 2 



dt 



(12) 



It appears that (M ee ) and (M t y P ) are time integrals of the return probability with various 
power-law weighting functions. Noting that P 7 (t) has the form Po(^)e _7i , all these quantities 



1 The lower bound of the time integrals is actually the mean collision time r e above which diffusion 
takes place. 



5 



can be written as integrals of a single function £(7, H) that we call the spectral function^ 
5( 7 , H) = J P °^ H K ~^dt = J^d^J C 7 (r, r, #)dr (13) 
The different magnetizations can be given in terms of the successive integrals of this function: 

(M„(B)) = -^±SV{>y,H) (14) 

Kv(H) = ^^p S{2 \^ H )\o (15) 
where 5 , W( 7 ) = L°° d 7 n..- XT d7iS(7i)- 



IV. EXAMPLES 

We first consider the case of a ring of perimeter L connected to one arm of length b, 
see fig.l. Such a geometry has been considered in the strictly ID case without disorder 
Buttiker 19"85| , |Buttiker 1994| , |ATtkermans et al. 1991| , [M"ello 1993|| . It is expected that since 



the electrons will spend some time in the arm where there are not sensitive to the flux, 
the magnetization will be decreased. From eqs. (0,0), the function C^(r,r,H) can be 
straightforwardly calculated on the arm and on the ring. Spatial integration gives: 

5(7, H) = — ln[- tanh — sinh - — h cosh cos47ry9] 

2 Lff, 

where ip = $/0o- $ is the flux through the ring, O = h/e is the flux quantum, jtd = 
(L/L^) 2 . (M ee ) and M typ are given by successive integrations over 7 according to eqs.([L4],|I^). 



They can be compared to the case of a single loop without arm. We find that they decrease 
when b increases and saturate to a finite value when b 3> L^. The m th harmonics of the flux 
dependence of (M ee ) is smaller by a factor (2/3) m and the m th harmonics of M typ is smaller 
by a factor (2/3) m / 2 . 

In the case of a ring connected to two diametrically opposite infinitely long arms, fig.l, the 
m th harmonics of (M ee ) is smaller by a factor (4/9) m and the reduction factor is (2/3) m for 



the typical magnetization. This result is relevant for the experiment of ref. [[Mailly et al. 1995 



where the magnetization is measured for open and connected rings. Moreover, in the limit 



2 This function S is related to the logarithm of the spectral determinant defined in ref. [An 



dreev et al. 1995 ] . The number variance for a closed system can be written as 

where s is the spin degeneracy and (3 = 1 in zero field, (3 = 2 if time reversal symmetry is broken. 
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b ^> L$, the magnetization should be unchanged if reservoirs are attached to the arms. We 
propose that single ring experiments with appropriately designed arms should be able to 
measure these reductions. 

We now turn to the case of an infinite square lattice whose average magnetization will 
be compared with the one of an array of isolated rings. The eigenvalues of the diffusion 
equation can be calculated for a rational flux per plaquette (p = Ha 2 /<f)Q = p/2q. a is the 
lattice parameter. Defining 77 = a/L^, we find that the spectral function £(7, H) is 

S(<y,E) = --]T((ln(4cosh7 / -^, / i))>> (16) 

1 8=1 

where (((. . .))) = J 27r |^ J 27r |^(. . .). £j(0, fi) are the solutions of the determinental equation 
detM = where the q x q matrix M is defined by M nn = 2 cos(4mry9 + 9/q) — e, M n>n +i = 
M n+ i 7n = 1 for n < q — 1 and M\ A = M* x = exp(z/i). This is the matrix associated 
to the Harper equation known to be also relevant for other related problems like tight- 
binding electrons in a magnetic field ||Hofstadter 1976fl or superconducting networks in a 



field [[Rammal et al. 1983|| . 

The expression ( |i~6D can be compared to the spectral function for a square ring of perime- 
ter L = 4a: 

5(7, H) = — ln(cosh 4r/ — cos 4inp) (17) 

From eq . (|i~4D , one calculates (M ee ): 

eD d q r°° 
(M^) = UPo^ttY, / {ln(4 cosh 77 -£;(#, /x))}7/dr/ 
n 2 q fr{ Jr, 

which is a continuous and derivable function of the field. It can be compared with the ring 
magnetization which can be cast in the form: 

, AeD f°° sin47T09 , 
(M ee ) = Up / — ^——Vdv (18) 

7T Jri COSh AT] — COS AlKfi 

The magnetization density is plotted on fig. 2 for the ring and the infinite lattice. It is 
seen that the network magnetization density is about 25 times smaller than the ring mag- 
netization. Considering that on the array of square rings already considered experimentally 



evy et al. 1990| , |Reulet et al. 1995|| , the distance between rings is equal to the size of the 
squares, the number of squares is four times larger when they are connected. One then ex- 
pects only a factor of order 6 between the magnetization of the array of disconnected rings 
and the lattice. 

We have also calculated the magnetization of a chain of rings connected with arms of 



similar length, an obvious generalization of the experiment done in ref. ||Mohanty et al. 1996 



and we find that when the rings are connected the average Hartree-Fock magnetization is 
reduced by a factor 3. 
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In conclusion, we showed that persistent current experiments do not necessarily require 
the rings to be disconnected and we suggest to perform magnetization experiments on arrays 
of connected rings. 



V. TRANSPORT 

We now come back to the calculation of transport properties in term of the spectral 
function S '(7), with application to the well-known case of an open wire as an example. From 
eq.(H), the dimensionless conductance g = G/(e 2 /h) is given by (for one spin direction): 



(5g) = -2 J P(t)- 



dt 
t d 

where td = L 2 /D is the diffusion time, inverse of the Thouless energy. This correction is 
easily written as a function of S^): 



it \ 9 1 d dS 
{5g) = 2 — — = 2— 

T D C/7 OX 



where x = jtd = (L/L 



*2 



Similarly the conductance fluctuation can be found directly as a function of the dif- 
fusion probability p(r,r',t) [|Argaman 1996|| . Then, using some properties of the diffusion 



propagator, it can be cast in the form, for one spin direction: 

tdt 



(Sg 2 ) = 12 J Pit)'- 



D 



or 

1 d 2 S _d 2 S 



(5g 2 ) = 12^— = 12- 



In the example of an open wire, the spectral function S is found to be, considering the 
appropriate boundary conditions: 

S(z) = -ln- Sinh ^ 



from which the weak-localization correction and the variance of the conductance fluctua- 
tions is immediately found. In particular, in the limit of complete phase coherence L v ^> L, 
i.e. x — > 0, the expansion of the function S(x): 



S{x) -> -- + 



2 

X X 



6 180 

immediately leads to the known universal values: 
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(Sg) = 2S'(0) = -1/3 

(Sg 2 ) = 12S"'(0) = 2/15 

The variance should be multiplied by two when there is time reversal symmetry. When 
is finite, the variance of the fluctuations is simply related to the mean weak-localization 
correction: 

(Sg 2 ) = -6-H-(6g) (19) 



In conclusion, we have written various thermodynamic and transport properties in term 
of a single function called spectral function, which is related to the spatial average of the 
return probability. We are thus able to calculate these physical properties for any type of 
network. 

Part of the work was done at the Institute of Theoretical Physics, UCSB, and supported 
by the NSF Grant No. PHY94-07194. 
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